Electronic interactions in fullerene spheres. 
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The electron-phonon and Coulomb interactions inCgO) and larger fullerene spheres 
are analyzed. The coupling between electrons and intramolecular vibrations give 
corrections ~ 1 — 10 meV to the electronic energies for Cqq, and scales as R~ 4 in 
larger molecules. The energies associated with electrostatic interactions are of order 
~ 1 — 4 eV, in Ceo and scale as R . Charged fullerenes show enhanced electron- 
phonon coupling, ~ 10 meV, which scales as R~ 2 . Finally, it is argued that non only 
■ CgQ, but also CgQ - are highly polarizable molecules. The polarizabilities scale as R 3 

and R 4 , respectively. The role of this large polar izability in mediating intermolecular 
interactions is also discussed. 
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I. INTRODUCTION. 



Fullerene molecules are being discovered in an atonishing variety of sizes and shapes. The 
most studied of them, C 60 shows a number of unexpected electronic properties, most notably 
the fact that, when combined with many metallic elements, it gives rise to superconducting 
compounds. Recently, a number a other spherical shaped carbon molecules have been found 
011 . In the present work, we discuss the most relevant intaractions which influence the 
highest occupied electronic states in these molecules. 

The analysis is based on earlier work 0,|J, (hereafter referred as I and II), which proposed 
a continuum scheme to analyze the properties of these molecules. The method makes use of 
the fact that simple, long wavelength approximations to the electronic spectrum of graphite 
planes can be made. We showed that it is possible, in this way, to write down analytical 
expressions for the highest occuppied, and lowest unoccuppied states of Ceo- Moreover, the 
method is more accurate for larger molecules. In I, a brief outline of its usefulness for the 
study of the coupling between electrons and lattice vibrations was also sketched. 

In section II, we describe the most relevant features of the electronic states of Cgo- The 
next section deals with the main properties of the coupling between electrons and phonons 
in neutral molecules. Section IV analyzes the electrostatic coupling between the electrons. 
Section V describes the new electron-phonon couplings which appear in charged molecules. 
In section VI we discuss the electronic polarizabilty of these molecules, and the interactions 
which may be relevant in a crystalline structure. The main conclusions are analyzed in 
section VII. 

II. THE ELECTRONIC STATES OF C 60 . 

The electronic properties of conjugated carbon compounds are determined by the states 
derived from the unsaturated 7r orbitals at each carbon atom. The very succesful Huckel 
model || uses this fact to describe a large variety of such materials, from benzene to conju- 
gated polymers, like polyacetylene, to graphite. The most important parameter in the model 
is the hybridization energy between orbitals placed at nearest neighbor carbons, which we 
will take as t ~ 2.2eV. Longer range hybridizations are neglected. The remaining three 
valence electrons in each C atom give rise to covalent o bonds. These bonds give rigidity to 
the lattice. Their spring constant, K, is also well known, which allows a simple approxima- 
tion to the vibration spectra of materials with only carbons, like graphite. We will assume 
for the highest frecuency, u opt = 3yjK/M = 0.20eV (M is the mass of a C atom). The 
electronic states derived from the n orbitals are affected by the lattice vibrations, because 
they modify the distance between nearest atoms. This coupling is described by the variation 
of the hybridization t with distance. We will take dt/dl = 4.5eVA . There is a greater 
controversy about the value of the on site Coulomb repulsion, when two electrons occupy 
the same n orbital. In isolated fullerene molecules, metallic serening is absent, and the long 
range electrostatic interactions dominate. We expect also no metallic screening at distances 
shorter then the intermolecule distance in Cqq crystals. In contrast, screening effects will be 
significant in onion like compounds, where the spheres are stacked inside each other. 

The Huckel model, when applied to a plane of graphite, gives rise to two electronic 
bands which intersect at two inequivalent points of the Brillouin zone. Near these points, 
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the spectrum can be written as: e k = hvp \k\, where vp — 3ta/h, and a is the distance 
between nearest carbon atoms. As discussed in I and II, this dispersion relation can be 
seen as arising from an effective long wavelength equation, which is nothing but the Dirac 
equation in (2+1) dimensions. The two inequivalent points of the Brillouin zone lead to 
two independent equations. The components of each spinor correspond to the amplitude of 
the wavefunction in each of the two sublattices in which the honeycomb struvture can be 
divided. 

This long wavelength approximation can be extended to spherical fullerenes. Naively, 
one would expect that the electronic states will be given by the solutions of the (2+1) Dirac 
equation on the surface of a sphere. The presence of pentagons in the lattice modifies this 
picture. When going around one of them, the two sublattices are interchanged. One of 
the components of the spinor associated with one point of the Brillouin zone is changed 
into the other component of the spinor associated with the second point. This effect can 
be included by means of an auxiliary gauge potential, whose magnetic flux threads each of 
the pentagons. Because of the mixing of the two electron flavors, this field is non Abelian. 
However, as discussed in II, a transformation is possible, which changes the two coupled 
Dirac equations, into two decoupled ones, which include an Abelian monopole of opposite 
charge, Q = ±3/2. 

These equations can be solved analytically, by means of suitably defined angular mo- 
mentum operators. The rwsulting spectrum has two degenerate triplets at e = 0, and a 
succession of higher multiplets, with degeneracies 2 x (2j + l),j = 2,3,... and energies 
€j = ±v F / R^J j(j + 1) — 2. R is the radius of the molecule. As extensively analyzed in I 
and II, this spectrum is a reasonably good approximation to spherical carbon fullerenes, 
becoming more accurate as the number of atoms in the molecule increases. In particular, 
for C 60 , the two degenerate triplets at e = mentioned above, can be identified with the 
lowest low lying molecular orbitals in the molecule. The fact that these two triplets are 
split in Ceo, and that they are not strictly at e = can be understood as arising from the 
next leading corrections to the Dirac equation. These corrections modify the spectrum of 
graphite away from the Fermi points. They contain additional derivatives. When included 
into our effective Dirac equation for the spherical fullerenes, they give rise to a coupling 
between the two triplets, and an overall shift. The fact that this interpretation is correct 
can be seen from the way in which these terms scale with the radius of the molecule. The 
splitting of the two triplets shows the expected R~ 2 dependence for all sizes, including Cgo- 
The behavior of the average energy is more involved, and the proposed scaling seems to set 
in only for large sizes (C 540 and beyond). 

In this way, we complete our description of the two lowest unoccupied triplets in C 60 
and related molecules. In lowest order, they correspond to the states at zero of the Dirac 
equation on the surface of a sphere, with an additional monopole of charge Q = ±3/2. Their 
wavefunctions are: 
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(1) 



The superscripts a and (3 refer to the two components of each spinor, and the subscripts 
label each of the six different wavef unctions. 

The residual interactions described earlier can be projected onto the subspace spanned 
by these six states. There is a global shift, which scales as k\t(Ro/R) 2 , and an hybridization 
between spinors of type a and b with the same angular momentum. It can be written as 
k 2 (Ro/R) 2 . Ro is the radius of C 60 , and the dimensionless coefficients k\ and k 2 can be 
deduced from the the spectrum of C 60 . Their values are 1.53 and 0.09 respectively. 

III. ELECTRON-PHONON INTERACTION IN NEUTRAL MOLECULES. 

We now describe the influence of the various lattice vibrations in the states given en eq. 
1. It is interesting to note first, that the hopping t, and any modulation of it, couple orbitals 
in different sublattices. To lowest order, the orbitals given by (1) represent electronic states 
which combine wavefunctions derived from one point in the Brillouin zone and one sublattice, 
and wavefunctions from the other point and the other sublattice. The hopping does not mix 
them, and that is why their energies are at zero. Thus, the only electron-phonon possible, 
for these triplets, arise from the modification by the phonons, of the next order terms, which 
were described before. This effect was not properly included in I. 

The states in equation (1) are delocalized thoroughout the sphere, and change little from 
site to site. Hence, they will only couple to long wavelength phonons. The graphite planes, 
which we use as starting point, have two accoustic and two optical branches near the center 
of the Brillouin zone, which we will consider separately. 

The accoustical modes are well described by the elastic theory of spherical shells, as 
mentioned in I. They can be parametrized in terms of fields which give the lattice displace- 
ments at each point in the sphere: u = (u r (9, <p),u e (9, 0), u^{6, <fi)). On general grounds, we 
can clasify them using vector spherical harmonics for the functions Uq and u^, and ordinary 
spherical harmonics for u r . For each value of I and m, the calculation of the eigenmodes is 
reduced to the diagonalization of a 3 x 3 matrix. The simplest vibration is the breathing 
mode, / = m = 0, u r = u, ug = = 0. This mode simply changes the radius of the sphere. 
We know that the lowest order coupling possible is through the modification of the electronic 
interactions within the triplets which decay like R~ 2 , as discussed in the previous section. 
Thus, its is straightforward to describe the coupling of the breathing mode to the lowest 
lying triplets given in eq. (1). To leading order in u/R, their mean energy and splitting 
depend on u as: 



e — h^u^ 
Ae = k 2 %u% 
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(2) 



For a general mode, the coupling to the electrons can only depend on the strain tensor, 
given by: 



ue,e = u r + d g u e 

U><t>,<f> "r T sin ^) T sin ^ ug 



(3) 



We now make the assumption that the way the strain tensor changes the electronic 
levels is through its only scalar contraction, ug t g + u^. Physically, this quantity describes 
the local variation in the area of the sphere. Generalizing the coupling to the breathing 
mode, discussed before, and particularizing in the changes in e, which is the largest effect, 
we write: 

dt f 

n el . phon[acc) = 2h— / {ue,e + (4) 

For an elastic vibration with a given I and to, ug t g + oc Y l m (0, <fi). Given the shape 
of the wavefunctions (eqn. 1), this result implies that only modes with I = 0, 1, 2 couple to 
them. 

We can now estimate the change in the electronic levels due to each lattice mode, using 
second order perturbation theory. From < u 2 >= h/(2Muj) (M is the mass of the molecule) 
the characteristic energy shift induced by a mode of frecuency u is: (k\dt / dl) 2 / \2Mu 2 ) . For 
Ceo, typical phonon frecuencies are in the range 500 - 1000 cm -1 . The associated shifts are 
of order 1-5 meV. The dimensionless electron-phonon coupling parameter, A, is < 0.1 per 
mode. 

In larger molecules, M ~ R 2 , and uo ~ it! -1 . Hence, the shifts will decay as R~ 4 , and 
A ~ i?~ 3 / 2 . Note that the number of accoustic modes which couple to the states studied 
here does not scale with the size of the molecule, because the characteristic dimension of 
the electronic wavefunction is always proportional to the radius of the sphere. 

The estimates discussed before point out to a weak coupling between the accoustical 
modes and the lowest lying electronic states of spherical fullerenes. This effect can be traced 
back to the low density of states of graphite near the Fermi level. For instance, the lowest 
order coupling between the triplets in eq (1) and phonons vanishes, but is finite for other 
multiplets at higher energies. Furthermore, the extended nature of these wavefunctions 
imply that few modes can couple to them. 

In molecules with one or two electrons in these states, a novel effect may arise. As 
discussed in I, the degeneracy of the allowed states gives rise to a Berry phase which mod- 
ifies the quantization rules for the lattice vibrations. The combined electron and phonon 
wavefunction changes its symmetry with respect to the case of no electron degeneracy. As 
a consequence, Raman and infrared modes are exchanged, an effect which can be verified 
experimentally. 

Finally, we discuss the coupling to optical modes. At each point in the sphere they can be 
characterized by its polarization, longitudinal or transverse, which can be defined as a two 
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dimensional vector field, n(#,0). We will neglect the out of the sphere modes. It is simple 
to calculate the coupling between this field and the spinors which describe the electronic 
states in planar graphite. It gives rise to a term in the hamiltonian: 

Hei-o P t = J ^anijj (5) 

where the cr's are the Pauli matrices. As in the case of the accoustical modes, this term 
does not modify, to lowest order, the triplets given in eq. (1). 

We can estimate the residual coupling by explicitly changing the hoppings in the discrete 
hamiltonian which describes Cqq, simulating the presence of an optical phonon. The sim- 
plest such mode alters the pent agon- hexagon and the hexagon-hexagon bonds in opposite 
directions. In order for the area of the sphere to stay constant, the pentagon-hexagon bond 
deformation should be close to one half that of the hexagon-hexagon bonds. The change in 
the average energy of the lowest triplets, with respect to such a modification of the hoppings 
is de/dt = 0.9. On general grounds we can expect this derivative to be of order unity. In our 
continuum formulation, such coupling can be incorporated into the formalism by a term: 

«~~ « - 2k 4 § / (*- + It + s®-) * r *- (6) 

where k% = 0.9. As in the case of accoustical phonons, this expression makes clear that 
only modes with small l,m values couple to the triplets considered here. 

The change in the energy of the electronic states is ~ (k^dt/dl) 2 /{2Mu 2 pt ). Taking 
^opt = 0.20eV, this energy is ~ 40meV. The dimensionless coupling A is ~ 0.2 — 0.3 per 
mode. These values for the energies and dimensionless couplings should be independent of 
the size of the molecule. Note that u opt does not scale with R, and M ~ R 2 . 

The combined coupling to all phonons compares well with more detailed calculations for 
^60 IHI- O ur definition of the dimensionless coupling A is dt/dl\J%/ (Mw 2 ) / (tuo), where uj 
is the frecuency of the mode under consideration. This is equivalent to the formulae used in 
H H|, in the limit when the coupling is local and the electronic bandwidth tends to zero. 



IV. COULOMB INTERACTIONS IN CHARGED FULLERENES. 

From equation (1) it is simple, although tedious, to calculate the different Coulomb 
integrals between orbitals. The needed integrals are discussed in the appendix. A variety of 
estimates are available in the literature, and our results come close to those reported in ref. 
@. In fact, if we neglect the splitting between the two triplets, we obtain, for Cg ~, exactly 
the same succession of multiplets as in that reference, in what is called also the continuum 
approximation. We believe that the scheme used in is equivalent to ours, in that the 
same wavefunctions are used. Thus, for Cg _ the lowest configurations are two 1 s + and 
singlets at 3e 2 /5i? 0) followed by two 3 p + and 3 p_ triplets, at 7e 2 /10i? - The splitting between 
these two sets of configurations is generated by the second order term which couples the one 
electron triplet states, discussed after eq. 1. An effective hamiltonian which describes each 
symmetry ( s and 3 p) is given in the appendix. These terms are also responsible for the 
splitting of the triplets in Cg . In the doubly ionized molecule, however, this term induces a 
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much weaker effect. The coupling of the two degenerate singlets is through an intermediate 
state at 8e 2 /5i? , an d the coupling between the triplets is mediated by another triplet at 
6e 2 /5R . In both cases, the splitting goes like Ae 2 /(e 2 /R ). For large fullerenes, Ae ~ R~ 2 , 
and this splitting decays like R~ 3 . Further details can be found in the appendix. 

Our method is sufficiently simple to allow us to calculate electrostatic energies of 
molecules with higher charge. For Cg , we find, at low energies, two 4 s multiplets. The 
electrostatic energy is 21e 2 /10i?o- The splitting between these two configurations is even 
smaller than in Cg ~, and goes like Ae 3 / (e 2 / Rq) 2 . For large fullerenes, this value scales as 
R~ A . In molecules with higher ionization, the splittings grow again. 

Finally, our simple form for the wavefunctions (1), allows us to characterize the charge 
distribution of the molecule. The two orbitals with angular momentum in (1) carry an 
electric dipole, eRo/2. It is easy to show that the two 3 p multiplets discussed for Cg ~ carry 
the same dipole. The existence of these large dipoles is an amusing consequence of the 
fact that the orbitals given in (1) simulate the state of electrons in a fictitous magnetic 
field, described by the monopole. These orbitals are not chiral invariant. A given angular 
momentum also induces an electric dipole. The terms responsible for the mixing of these 
orbitals restore the chiral symmetry. As mentioned before, the influence of these terms is 
significantly reduced in Cg~ ~ and Cg , and in larger fullerenes. The relevance of these 
dipoles for intermolecular interactions will be discussed in section VI. 

The scale of electrostatic energies, given the radius of C 60 , -Ro ~ 3.5A is large, ~ 3 
eV. We have not investigated to which extent the level structure of the spheres needs to 
be modified. For larger fullerenes, both the separation between one electron levels and the 
Coulomb energies scale in the same way, ~ R~ x , so both remain comparable for all sizes. 

An alternative method, used to classify the multiplets of the C 60 molecule is to use 
spherical harmonics and crystal field terms pH|jn[ ]. This scheme is equivalent to use a 
continuum approximation for the bottom of the graphite band (the point k = in the 
Brillouin zone), and to fill all the levels up to the required charge state. Our method has 
the right multiplicity and symmetry for the relevant orbitals near the neutral state, as it 
starts from an approximation which describes graphite near its Fermi energy. Moreover, the 
separation between levels, which in other schemes arises from crystal field effects, appears in 
a natural way, and its scaling as function of size can be easily investigated. In that respect, 
we differ somewhat from the assumptions about scaling in ref. |T . 



In onion shaped materials, the energy levels are affected by the tunneling between con- 



tiguous layers. The energy of these processes is, in graphite, ~ 0.3 eV [12]. The Coulomb 



interactions will be more drastically changed, by the appearance of screening effects. 



V. ELECTRON-PHONON INTERACTION IN CHARGED FULLERENES. 

In charged fullerenes, the deformations of the lattice induce rearrangements in the charge 
distribution, which modify the Coulomb energies. As discussed in the previous section, 
these energies are important, so that this effect needs to be taken into account. This effect 
is unique to systems like the doped fullerenes, which combine localized electronic states 
and inhomogeneous charge distributions. Hence, other conjugated compounds, like graphite 
intercalates, offer no clue to the strength of this coupling. 
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These modifications in the charges over scales comparable to the size of the molecule can 
only be induced by long wavelength accoustical phonons. Optical modes redistribute the 
charge within each unit cell, but do not alter the overall charge distribution. 

The leading effect of a deformation of the lattice is to change locally the area element. 
The charge density in that region is then modified. Again, we can discuss this effect taking 
the simplest lattice vibration, the breathing mode as an example. The area of the sphere 
is expended or contracted, and the Coulomb energies change due to the variation in the 
radius of the sphere. There is no effect in Cg . In Cgo _ , the change of energy associated to 
a vibration of amplitude u is ~ e 2 u/R 2 r For a general mode, the charge density at a given 
place decays as the area element increases. The local expansion induced by a given mode 
is ue,e + u^. Thus, we have to insert this expression in the Coulomb integrals between 
different orbitals. As in the case of the standard coupling to accoustical phonons, this result 
imply that only vibrations with / < 4 can couple to the triplets given by eq. 1. 

We can make an estimate of the energy shift in the electronic levels in the same way 
as in the case of the modification in the hybridization due to the vibration. The only 
difference is that the factor {k\dt/dl) 2 has to be replaced by (e 2 / R^) 2 ■ The second factor 
is larger by a factor of 8 - 10. Thus, the energy scales at wich this coupling can play a 
role are ~ 20 — 40meV. The values of the dimensionless constant A will become ~ 0.8 — 1. 
These interactions are significantly larger than the standard electron-phonon couplings. The 
energies should scale as R~ 4 for large molecules, and the values of A as R~ 3 / 2 . 

VI. INTERACTIONS BETWEEN DIFFERENT FULLERENE MOLECULES. 

So far, we have discussed the various electronic interactions within a given fullerene 
molecule. It is interesting to analyze which ones may play a role when different molecules 
are close to each other. As the most important influence that we have found is due to the 
Coulomb interaction, we will analyze further its role in this case. 

We will assume that charge that the fullerenes may have is compensated by some neu- 
tralizing background, so that the only residual interaction is due to the Van der Waals forces. 
Then, the relevant magnitude is the electrical polarizability of the molecules. As mentioned 
in section IV, Cg is highly polarizable, a result also discussed in ref. 0. Also, Cgo~ has a 
large polarizability, if the molecule is in a 3 p state. An estimate of the polarizabilities of 
various charged states of fullerenes is given in table I. These polarizabilities are due to the 
extra electrons which are present in charged Cgo molecules. The polarization of the closed 

9 3 

shell in neutral Ceo has been estimated to be ~ 80A ~ 300 a. u., that is, significantly 
smaller than the values we find for Cg and Cg^ - . 

The high polarizability of Cg ~ can play a role, if the interaction between neighboring 
molecules overcomes the energy difference betwwen the 1 s and the 3 p configurations. This 
difference is e 2 /(10i?o) ~ 0.5eV, as discussed in section IV. The van der Waals interaction 
between two molecules at distance D is: ~ e 4 i?o/(4-D 6 e p ), where e p is the splitting between 
the two 3 p configurations, ~ 2Ae 2 /(e 2 /Ro) (see the appendix for details of the calculation). 
Hence, the van der Waals energy is ~ 5e 6 -Rg/(4/} 6 Ae 2 ), where Ae ~ 0.40eV, as discussed 
in section 2. Even for D = 4i? , this energy is ~ 0.2eV, and increases rapidly at shorter 
distances. Thus, it is not unconceivable that Cg ~ will exist in its 3 p state in a crystal, and 
neighboring molecules will interact strongly through their mutual dipoles. 
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VII. CONCLUSIONS. 



We have investigated the various interactions which may influence the electronic prop- 
erties of a fullerene molecule. The relevant energies, and the scaling as function of molecule 
size are given in table II. 

Our results suggest that the most relevant effects will arise from the Coulomb interac- 
tions. The weakness of the electron phonon coupling can be related to a similar effect in 
graphite, due to the vanishing density of states near the Fermi points. Phonons can only 
play a significant role in charged molecules, where they modify the charge distribution. 

On the other hand, long range Coulomb interactions have no counterpart in other con- 
jugated compounds. The most studied charged (through doping) carbon systems, graphite 
intercalation compounds and poliacetylene, do not show these features. In both cases, the 
relevant electronic states are delocalized, and are not susceptible to charging effects. 

We have pursued further the role of Coulomb interactions in crystalline systems. The 
most striking effect that we have found is the extremely high polarizability of Cg _ molecules. 
This property leads to significant intermolecular interactions at distances ~ 3 — 4 times the 
radius of the molecule. The scale of these couplings ~ 0.3 — 3eV, can easily overcome most 
other efffects in a doped Ceo crystal, and may be even responsible for the unusual dipolar 
moments found in neutral C 60 P^| . Its role in other properties, like superconductivity, is 
being further investigated. 

CALCULATION OF COULOMB INTEGRALS. 

The Coulomb integrals required to evaluate the energies of the various configurations 
of charged Cgo can be obtained from the wavefunctions given in eq. (1). We denote the 
six orbitals in that equation as: \a, +1 >,|a,0 >, |a, — 1 >,|6, +1 >, |6, > and \b, — 1 >, 
where we are omitting spin indexes. The Coulomb potential can be expanded in spherical 
harmonics, l/\fi — r 2 | = 4:iiJ2i 1/(2/ + 1) Em Y l ~ m (ri)Y l m (r2), where we are using the fact 
that |fi| = |r 2 | = 1, and we use units such that R Q — 1. Then, the seven different Coulomb 
integrals are (in units of e 2 /R ): 

V a +l,a+l;a+l,a+l = ^a+l,a+l;6-l,6-l 

^a+l,a+l;a0,a0 = ^a+l,a+l;bO,feO 
v a+l,a+l;a-l,a-l = v a+l,a+l;b+l,b+l 
v a+l,aO;a+l,aO 
V a +l,aO;b+l,bO 
Va+l,a-l;a+l,a-l 
VaO,aO;aO,aO = V a 0,a0;60,ft0 

The remaining integrals can be obtained by exchanging the a and b indices. 
As mentioned in section IV, this procedure leads to the multiplet energies for C^ _ dis- 
cussed in ||. The wavefunctions with 1 s symmetry are: 

|1 > = ^(|o + l |,a- 1 |> -|a0 |,a0 J> +|a - 1 t, a + 1 |> ) 
|2> = X (|a + l T,&-1 I> -|a + l t> +|a-l T,&+1 I> 



1 

_ 2 ii 

„50 



23 



(7) 
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- \a - lj, b + 1 t> ~\aO t, bO |> +|aO |, 60 |> 
3 > = ^(|6+1 t,6- U> -|60 T,W |> +|6- 1 T,6 + U> ) 



(8) 



These three states are mixed by the term in the Cqq hamiltonian which gives rise to the 
splitting of the two triplets, Ae. Thus, we end up with the matrix: 



Hi. = 



/ 3ei Ae n \ 

/ 5Ro V? \ 

Ae 8ef_ Ae 

in A£ 3ei / 

V U V2 5R J 



Similarly, the states with 3 p symmetry are: 
U > 



= \a+ 1 T,a0 |> 
2> = ^(|a + l t,60 T> +I&+1 T,o0t> 
3> =|6 + 1T,60T> 



and the corresponding effective hamiltonian is: 



Ha 



' lORo V2 U 

Ae 6ei Ae 

V2 5R y/2 

in ^ ^- I 

\ y/2 lORo / 



(9) 



(10) 



(11) 



For the purposes of the discussion in section VI, it is important to realize that the states 
|1 > and 1 3 > in (10) carry an electric dipole, ±eR /2. 
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TABLES 

TABLE I. Typical energy scales associated with the various interactions discussed in the text, 
for Ceo, and scaling of each of them as function of molecule radius, for larger fullerenes. 



Interaction 


Energy 


Scaling 


Electron-phonon (accoustic) 


0.004 eV 


ir 4 


Electron-phonon (optical) 


0.04 eV 


R° 


Coulomb 


3 - 6 eV 


R- 1 


Electron-phonon (in charged molecules) 


0.1 eV 


R-* 



TABLE II. Polarizability (of the valence electrons) in different charge states of Cgo (in atomic 
units), and scaling behavior for larger molecules. 



Charge state 


Polarizability (in a. u.) 


Scaling 


Neutral ( 1 s ) 





R° 


le"( 2 p) 


10 3 


i? 4 


2 e" ( x s ) 





R° 


2 e" ( 3 p ) 


10 5 


R 5 


3 e" ( 4 s ) 





R° 
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